Local power fluctuations in two-dimensional turbulence. 
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The statistics of power fluctuations are studied in simulations of two-dimensional turbulence in 
both inverse (energy) and direct (enstrophy) cascade regimes from both Lagrangian and Eulerian 
perspectives. The probability density function (PDF) of the appropriately defined dimensionless 
power is strongly non-gaussian with asymmetric exponential tails. This distribution can be mod- 
eled by the distribution of the product of correlated normal variables allowing a derivation of the 
asymptotics of the tails. The PDF of the dimensionless power is shown to exhibit an empirical Fluc- 
tuation Relation. An expression for the entropy production rate is deduced from the asymptotic 
form of the power PDF and is found to agree very well with the measured entropy rate. 

PACS numbers: 47.27.Gs 



INTRODUCTION 

Thermodynamics provides a theoretical framework for 
describing the macroscopic properties of equilibrium sys- 
tems which are largely independent of microscopic de- 
tails. Much research in recent years has searched for 
universal features in non-equilibrium systems to help de- 
velop macroscopic theories of non-equilibrium statistical 
physics that are similarly insensitive to microscopic de- 
tails. One important branch of this research has focused 
on the Gallavotti-Cohen Fluctuation Theorem [l]] and 
related results 0, H, 01 which are collectively referred to 
as Fluctuation Relations. A Fluctuation Relation (FR) 
is formulated for the PDF of some quantity, Xr, de- 
rived from the entropy production or energy dissipation 
in a system driven far from equilibrium. Xr is obtained 
by averaging a physical quantity, x{t), typically the en- 
tropy produced or energy dissipated along a phase space 
trajectory of the system, over a time interval [t,t + r]: 
Xr = /j x{t') dt' . Xr is positive on average but, 
may fluctuate sufficiently about its mean that negative 
fluctuations are observable. An FR quantifies the rela- 
tive probability of observing a negative fluctuation over a 
given time interval compared to the probability of observ- 
ing a positive fluctuation of the same magnitude. The 
ratio of probabilities takes the form: 



n{-Xr) 



(1) 



where E is a constant, independent of the averaging in- 
terval, T. In the original formulation of x{t) is the 
entropy produced along a trajectory in phase space of a 
microscopically reversible chaotic system satisfying a cri- 
terion called the "Chaotic Hypothesis" . In [ij , Eq. ([J) is 
a theorem in which E = 1 and there are no adjustable 
parameters. Attempts to probe the validity of these hy- 
potheses by experiments and numerical simulations have 



struggled with the fact that it is difficult to find systems 
for which the assumptions of the theorem, are under con- 
trol. See 5] for one such experimental test and further 
discussion. 

Notwithstanding this difficulty, several experimental 
and numerical studies have studied the analogue of Eq. ([T]) 
in situations where the result may not be expected to 
hold a-priori. See Q and the references therein. Despite 
the fact that, in almost all cases, these systems were dis- 
sipative and follow microscopically irreversible dynamics, 
they have succeeded in yielding an empirical FR. Here, 
empirical means that the parameter E, which is typically 
not 1, cannot be deduced a-priori from known properties 
of the system (with the possible exception of the granular 
experiment discussed in [7] ) . Not all of these experiments 
have provided equally good tests of the FR. In some ex- 
periments the fluctuations were Gaussian, in which case 
Eq.([T]) is immediately upheld. While it may not be triv- 
ial to prove rigorously the Gaussianity of fluctuations in 
an interacting system, Eq.([T]) is considered to be most 
interesting when the PDF of Xr is non-Gaussian. Other 
experiments could not access sufficiently large negative 
fluctuations to provide a statistically convincing test. Ef- 
forts are ongoing to try to find good systems on which to 
explore the seeming ubiquity of the FR. 

From the theoretical perspective, it has been pointed 
out that any Langevin system exhibits a fluctuation re- 
lation for a properly defined entropy production (al- 
though such a definition need not be unique) . While one 
could argue that some of the aforementioned experiments 
could be sensibly modeled by Langevin dynamics, most 
of them measure work or power rather than entropy pro- 
duction and the relation between the two is not always 
obvious. Indeed it has been shown analytically 
that power injection fluctuations in driven Langevin sys- 
tems generally do not exhibit a Fluctuation Relation. It 
has also been suggested ll| that the FR follows as a 
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property of large deviations having generic Kramer func- 
tions, that is, those which are differentiable near zero. 

In this letter, we contribute to this ongoing discussion 
in two ways. Firstly, we present results for a new system 
which exhibits an empirical FR - local power fluctuations 
in turbulence. While it may not help much to simply add 
another system to the list discussed above, we argue that 
this system offers considerable advantages as a model 
system compared with some of the others. It is easy 
to realise experimentally and numerically, has strongly 
non-gaussian fluctuations and exhibits abundant nega- 
tive fluctuations. Furthermore, we deduce an expression 
for the entropy production rate from known properties 
of the power distribution that agrees very well with the 
empirical results. 

Secondly, we provide an analytic model of the PDF 
of the power fluctuations themselves (before time aver- 
aging) based on the PDF products of normal variables. 
This PDF and its asymptotic properties match very well 
with our empirical PDFs and will find many applications 
beyond the present study. It is worth remembering amid 
the current flurry of interest in FRs, that the PDF of the 
injected power itself is probably of more direct physical 
and practical interest than the PDF of its time-averaged 
counterpart. On this basis, the results we present on the 
structure of the power PDF itself are of more general 
interest, beyond the topic of FRs. The layout of the pa- 
per is as follows. We first provide a brief discussion of 
two-dimensional turbulence and explain the power mea- 
surements in the Eulerian and Lagrangian frames. We 
then describe the structure and asymptotics of the PDFs 
of the power and compare with the empirical PDFs ob- 
tained from our data. Finally we present data demon- 
strating an empirical FR and explain the original of the 
value obtained for E. 



TWO-DIMENSIONAL TURBULENCE 

The considerable differences between the physics of 
turbulence in 2D compared to 3D can be traced to the 
existence of an additional inviscid invariant (in addition 
to the kinetic energy), known as enstrophy. As a result 
of the dual conservation of energy and enstrophy in 2D, 
there are two distinct scaling regimes in 2D turbulence 
in the inviscid limit. If energy and enstrophy are in- 
jected into the flow at a given scale. If, the enstrophy 



tends to be transferred to scales smaller than 



This 



process is refered to as a direct cascade. The energy, on 
the other hand, tends to be transferred to scales larger 
than If. This phenomenon, referred to as an inverse cas- 
cade is specific to 2D. According to the classical theory 
of Kraichnan [l^ . the energy spectrum in the direct 
cascade regime scales as where k is the modulus of 
the Fourier space wave-vector whereas the energy spec- 
trum in the inverse cascade regime scales as k~^^^. For 



detailed review of the theory and phenomenology of 2D 
turbulence see [ist and the references therein. 

Although we have studied the statistics of the en- 
ergy injection rate in both the direct and inverse cascade 
regimes, in this letter we primarily discuss results for the 
inverse cascade regime. The analysis performed for the 
direct cascade regime was observed to be essentially sim- 
ilar to the inverse cascade results. The only related study 
which we are aware of is lJ| where the statistics of the 
global (i.e. averaged over the entire system) power have 
been studied experimentally in 3D. Here, we study the 
local power in 2D as we now explain. The instantaneous 
rate of energy input at a point, x in the fluid is given by 
the scalar product of the instantaneous fluid velocity at 
X, v(x, t), and the instantaneous body force, f(x, t), act- 
ing on the fluid at x. To construct the PDF of -f there are 
two approaches. One can fix the time and measure the 
V and f at different points and construct the PDF over 
the values of the local power at different spatial points. 
This is called an Eulerian measurement. Alternatively 
we may follow the trajectory, x(t), of a particular fluid 
element as it is advected by the turbulence and measure 
the force, f(t) = f(x(t),t), velocity, v(t) = v(x(t),t) 
and power p{t) = v(t) • f (t) experienced by this element. 
We then construct the PDF over the values of the power 
at different times. This is called a Lagrangian measure- 
ment. We have measured the power distribution in both 
the Eulerian and Lagrangian frames but focus more on 
the latter. The variance of the velocity is slighly larger in 
the Lagrangian frame but once this has been taken into 
account, the power statistics look very similar for both 
cases. For a careful discussion of the similarities and dif- 
ferences between the Eulerian and Lagrangian frames in 
2D turbulence see 
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The forcing used in our simulations, models that used 
to excite turbulence in electromagnetically driven fluid 
layers, a popular method of generating 2D turbulence 
in the laboratory [3]. We simulated a biperiodic do 



main of size 27r. For details of the simulations see [17 1 . 
Aside from the issue of experimental relevance, the rea- 
son to choose the electromagnetic forcing protocol over 
the stochastic forcing protocols traditional for numerical 
simulations of turbulence is that we do not wish to in- 
troduce an additional stochasticity into the system. All 
of the power fluctuations should come from the intrin- 
sic stochasticity of the flow. To generate a direct cas- 
cade, we forced at large scales using a distorted magnetic 
fleld perpendicular to the layer with characteristic wave 
number 5. To generate an inverse cascade, we forced at 
smaller scales with a magnetic field having characteris- 
tic wavenumber 35. For further details of the turbulence 
diagnostics and spectra see 17[. The fluid was driven 



with a direct current in the a;-direction. The only compo- 
nent of the force is then in the y direction so that p(t) is 
just given by a simple product of y-components of force 
and velocity rather than a scalar product. Results for 
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(A) Force PDFs 



(B) Velocity PDFs 
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FIG. 1: PDFs of the power normalized by the product of 
standard deviations of force (ct/) and velocity {a^) for both 
the Inverse (circles) and Direct cascade (asterisks) regimes in 
the Lagrangian frame as well as the Inverse Cascade regime 
in the Eulerian frame (squares) . The solid line is the asymp- 
totic behavior predicted by Craig's XY-model in Eq. ^ for 
the power distribution. The insets show the corresponding 
(A) force PDFs and (B) velocity PDFs for the Inverse and 
Direct cascade regimes in Lagrangian framKF2007e, the In- 
verse cascade regime in Eulerian frame and their respective 
Gaussian fits. Please note the force and velocity PDFs have 
been normalized by their respective standard deviations. 



alternating current are discussed in One should ap- 

preciate that although the forcing is deterministic (the 
disorder of the underlying magnetic field is quenched), 
in the Lagrangian frame, it is sampled by a random tra- 
jectory so that the force acting on a Lagrangian particle 
is effectively stochastic. Similarly, from an Eulerian per- 
spective, the force is a random variable in the sense that 
it takes a distribution of values at different spatial points. 

Fig. [1] shows plots of the PDFs of the dimcnsionless 
power for different regimes. The non-dimensionalisation 
has been done by normalising with the standard devi- 
ations of the force and velocity, P ~ p/{ay(Tf). This 
normalisation quantifies our earlier assurances that the 
choice of inverse versus direct cascade or Eulerian ver- 
sus Lagrangian frame is unimportant. The shape of the 
power distribution is practically the same once the data 
have been rescaled by these standard deviations. We 
should point out at this point that ver y s imilar PDFs 
have been observed in similar contexts [lil . flil [20| . It 
turns out that the Lagrangian force and velocity are 
quite close to Gaussian in all cases (see also [2l|). The 
PDFs, again rescaled by the appropriate standard devi- 
ations, are shown in the insets for Fig. [1] along with a 
Gaussian indicator curve. The measured standard devi- 
ations and instantaneous correlation coefhcient between 
v{t) and f{t) for the inverse and direct cascades are given 



TABLE I: Statistical parameters of the Lagrangian force and 
velocity. 



Simulation 


^} 




P 


Inverse cascade 
Direct cascade 


0.082 
0.338 


0.047 
0.589 


0.110 
0.144 



in Table U 



LOCAL POWER AND THE FLUCTUATION 
RELATION 



The power PDF, is strongly non-gaussian as can be 
seen from Fig. [1] For both the direct and inverse cas- 
cades, it is strongly cusped at with asymmetric ex- 
ponential tails. Given that p{t) = Vy{t) fy{t) and both 
Vy (i) and fy (t) are close to Gaussian, one may ask if the 
non-trivial features of the power distribution can be un- 
derstood by considering the distribution of the product 
of two correlated normal variables [l3| [l^. The dis- 
tribution of the product, P = vf, of two joint normally 
distributed variables, v and /, with zero mean, standard 
deviations c/ and ay and correlation coefficient p was 
considered by Craig The PDF of the product is: 



nxy(P) 



VA+A- 
2^ 



'dw 



(2) 



where A* — {axcry{l ± p))~ ■ In general, this integral 
cannot be further simplified. The asymptotic behaviour, 
however, is relatively easily calculated. We find that the 
tails are indeed almost exponential with the positive tail 
decaying more slowly than the negative tail. Full details 
are given in [l7| : 



A+A- c--^^^ 
^(A++A-) VP 



A+A- 
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/|P| 



p > 
p < 



(3) 



FigjT] super-imposes this asymptotic behaviour on the 
empirical PDFs for the measured values of cr^, Uy and 
p for each cascade regime tabulated in Table |T1 Note 
that we plot the normalized power P = p/ {a^af) for the 
asymptotic behavior predicted by Eq. ([3]). The agree- 
ment is excellent. This agreement is a result of the fact 
that the values of p are close for these situations. Other 
values of p for other systems may lead to a greater or 
lesser degree of asymmetry. Near zero, where the asymp- 
totic expression in Eq. ([3]) is not valid, Eq.(l2]) can be 



17|. Although our 



shown to be logarithmically singular 
data is insufhcient to resolve such a weak divergence, this 
explains the cusp-like structure of the empirical PDFs 
near zero. We conclude from this analysis that most of 



4 



c 




Q_ 

c 

c 
Q_ 




FIG. 2: PDFs of the coarse-grained dimensionless power (Pn) 
as defined in Eq. ((S} for the Inverse cascade regime. The 
PDFs are plotted for the dimensionless averaging time n = 
t/tc = 0, 5, 10, 15, and 20 respectively. The coarse-graining 
leads to narrowing of the PDFs with increasing n as expected. 
However they retain the non-gaussian features of the non- 
coarsegrained distribution. 



the interesting features of the power PDF can be under- 
stood in statistical terms. Dynamical information enters 
only through the values of ct/, ay and p which, we should 
emphasise, are measured rather than calculated from first 
principles. 

We note that Eq.Q shows that the tails of the power 
PDF satisfy a fluctuation relation of the form Eq.(IT]). 
Taking the ratio of the right to the left tail from Eq. ([3]) 
results in an expression for S: 



2p 



(4) 



The entropy rate S is now expressed purely in terms 
of the correlation coefficient p , which is known to decide 
the degree of the asymmetry of the power distribution. 



see [17[ . It is natural to expect that the time integrated 
power should satisfy a similar relation with a comparable 
value of S. 

We test the FR for the Lagrangian power in the in- 
verse cascade regime since this is the regime for which 
we have the most data. We coarsegrained the measured 
power over time windows which are longer than the cor- 
relation time, Tc, of the Lagrangian power signal. These 
correlation times are quite short since Lagrangian quan- 
tities tend to decorrelate very quickly. For the inverse 
cascade we measured Tc = 0.144 and for the direct cas- 
cade, Tc = 0.084. These compare to large eddy turnover 
times of about 22 and 10 respectively. We then define 
the non-dimensionalised coarse-grained power over a non- 



FIG. 3: Plots for the Left Hand Side of Eq. © for differ- 
ent coarse-graining intervals n = 0, 5, 10, 15, and 20. The 
solid line is the theoretical prediction for E deduced directly 
by applying the XY-model to the non-coarsegrained power 
distribution. As is evident, the theoretical prediction for E 
is in excellent agreement with the observed empirical Fluctu- 
ation Relation. The inset shows the functional dependence 
of the entropy rate E to the correlation coefficient p which 
essentially decides the asymmetry of the power distribution. 



dimensional time-interval of length n = t/tc'. 



Pn{t) 



1 1 



t + T 



P{t')dt'. 



(5) 



Here, t' is a dimensionless time normalised by the cor- 
relation time Tc- In fluctuation relation analyses, it is 
customary to normalize the signal of interest by its long- 
time average. However, in the present case, we normalise 
the power p by Pq = a^af. From Eq. [2] we see that this 
is an attractive normalisation since it reduces the under- 
lying PDF to a function of p only. The statistics of Pn 
are stationary in t since the underlying turbulence had 
reached a stationary state before we started to gather 
power statistics. The coarse-grained PDFs are shown 
in Fig. [2] for dimensionless averaging times ranging from 
n — t/tc = to 20. They clearly narrow as one inte- 
grates out more and more of the large fluctuations but 
retain their non-gaussian character as the averaging pro- 
ceeds. If the FR is satisfied, then Eq([T|) suggests that 
all the data should collapse according to the following 
formula: 



-log 
n 



n(+p„) 



H(-p„ 



(6) 



Additionally, if the asymmetry of the power distribu- 
tion observed for the non-coarsegrained power PDF is 
retained after coarsegraining, one expects the left hand 
side of eq. ^ should scale with the slope predicted by 
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eq. (13]). This collapse of the data is demonstrated in 
Fig [3] for the Inverse cascade. The value of S obtained 
from Eq. ^ after substituting the appropriate value for 
p from Table U is 0.22. The corresponding value of S 
obtained from the empirical FR plot in Fig. [3] is 0.217 
for the inverse cascade regime. The solid line in Fig [3] 
is the theoretical prediction following eq. [4] and is in 
excellent agreement with the empirical Fluctuation Re- 
lation. Note that owing to the coarsegrained power and 
time being non-dimensional, the entropy rate S, is also a 
dimensionless quantity. In the inset to Fig. [3] we plot the 
expected functional dependence of E with respect to the 
correlation coefficient p. The value of p and hence the 
asymmetry of the power PDF changes with the strength 
of the turbulence, which is normally quantified in terms 
of the dimensionless Reynolds number (Re). Depending 
upon how p varies with respect to i?e, one is then in a 
position to understand how the entropy rate S varies as 
a function of the Reynolds number. Such an analysis is 
however beyond the scope of the present work. 



SUMMARY 

In conclusion, we have demonstrated an empirical 
FR for the statistics of Lagrangian power in turbu- 
lence. It is an attractive system in which to study 
non-equilibrium fluctuations since the fluctuations are 
strongly non-Gaussian with an abundance of strong neg- 
ative fluctuations. In addition, we have phenomenolog- 
ically established the analytic form of the underlying 
power PDF. Its qualitative features can be understood in 
terms of the statistics of products of normally distributed 
variables. Most importantly, we are able to arrive at a 
theoretical expression for S that is in excellent agreement 
with the empirical results. 
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